Fragmented many-body ground states for scalar bosons in a single trap 
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We investigate whether the many-body ground states of bosons in a generalized two-mode model 
with localized inhomogeneous single-particle orbitals and anisotropic long-range interactions (e.g. 
dipole-dipole interactions), are coherent or fragmented. It is demonstrated that fragmentation can 
take place in a single trap for positive values of the interaction couplings, implying that the system 
is potentially stable. Furthermore, the degree of fragmentation is shown to be insensitive to small 
perturbations on the single-particle level. 



Fragmentation in many-body states of interacting 
bosons is denned by the single-particle density matrix 
having more than one macroscopic eigenvalue 0,0], lead- 
ing to macroscopic occupation of more than one state. In 
the continuum, and for repulsive contact interactions, the 
Fock term in the total energy favors a single condensate; 
a negative contact interaction implies the collapse of the 
gas before ground-state fragmentation sets in [3j . On the 
other hand, in spatially well-separated systems, fragmen- 
tation can be obtained for deep double wells [4j, or in its 
periodic extension Mott state in optical lattices [al. Fur- 
ther examples for fragmentation have been found [6|], e.g., 
when internal degrees of freedom and sufficient symmetry 
in the interaction are present Q , in the Richardson pair- 
ing model @, in the crossover to fermionization and for 
excited states jjj, and in rapidly rotating gases 1^, 12 1- 
However, no ground-state fragmentation of a scalar Bose 
gas in a single trap has been reported. 



So far, only contact interactions between the bosons 
and quasi-homogeneous or spatially periodic systems 
were considered. The question then arises whether for 
general inhomogenity, nonlocality and anisotropy, i.e. 
given the localization in a single trap with inhomoge- 
neously distributed single-particle orbitals (modes) and 
for partly positive and negative long-range interactions 
(of which the archetype is the dipole-dipole interaction) , 
ground-state fragmentation can occur. In the following, 
we consider for simplicity and easy comparison with the 
double-well case, that only two single-particle orbitals 
are occupied Examples are gases strongly confined 
in a given direction, e.g., close to the quasi-2D limit; the 
question then arises whether fragmentation occurs with 
respect to that strongly confining direction. A two-mode 
approximation has, furthermore, the important benefit 
that the many-body states of the system can be found 
essentially exactly numerically (and for special cases an- 
alytically). We demonstrate that by tuning the four in- 
teraction matrix elements relative to each other, various 
many-body states - coherent and fragmented states, as 
well as coherent superpositions of degenerate macroscop- 
ically distinct quantum states [l3| - can be accessed in 
a single trap. In addition, the fragmented states are not 
susceptible to decay to a non-fragmented (coherent) state 



because of a perturbation on single-particle level 

nana. 

due to the fragmentation being based to the values of the 
interaction couplings. 

We begin with a general quadratic plus quartic Hamil- 
tonian for two interacting modes, 



H 



e oja + eiajoi 



"2 



h.c. 



2 
A 



.4 



— a[a[aiai 



(1) 



A3 /„ t ~t 
~2 {%% 



a\ai 



— a{aia [ a . 



The interaction coefficients are given by A\ = 
Vbooo,A 2 = Vim, A3 = Vuoo = Vbon (taken to be 
real), and A 4 = Vbioi + Vioio + V1001 + ^0110 , where 
V m = Jd 3 rJ d 3 r'^*(r)^*(r')Vi nt (r-r')«' fe (r')^(r-). 
It is of major importance for the discussion to follow that 
we include pair-exchange between the two modes due to 
scattering of pairs of bosons oc A3, in addition to the 
standard density-density type terms oc A\, A-i, A^. We 
stress that this term is absent in the homogeneous con- 
tinuum due to momentum conservation. We also note 
that the independence of the coefficients Aj makes (U) 
different from rotationally invariant interaction Hamilto- 
nians in gases with internal degrees of freedom, where 
the modes are components of a spinor 0, 0|. For Fourier 
transformable interaction potentials we have, by convolu- 
tion, Vijki = (2^W / d 3 kpu(-k)V in t(k)pjk(k), where the 
Pjk(k) are Fourier transforms of \E , *(r)^f , / £ (r). The single- 
particle energies a = J d 3 r[^\V^i\ 2 + Vt rap |*i| 2 ] are 
given by kinetic plus trap contributions. Finally, the 
Josephson-type coupling fl connects the single-particle 



states and corresponds to a tunneling rate 

We emphasize that differing magnitudes and signs of 
the A can stem both from the interaction and the form 
of the single particle orbitals. This only simplifies in 
the continuum limit (the orbitals are plane waves then), 
where the interaction couplings directly reflect anisotropy 
and long-range character of the interaction. Note that 
for contact interactions, all Aj negative necessarily im- 
plies collapse for large systems, occurring by the fact 
that choosing ever more singular single particle orbitals 
is energetically advantageous. On the other hand, by ad- 
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justing the single-particle orbitals (the trapping) in con- 
junction with in general anisotropic interactions like the 
dipole-dipole interaction, the magnitude and sign of the 
interaction coefficients Aj can be engineered [2(j j. The 
Hamiltonian |T]) therefore represents a minimal model 
to investigate whether ground-state fragmentation takes 
place for potentials and traps with arbitrary anisotropy 
and shape, leading to essentially independent interaction 
couplings Ai. In a single trap, the two terms in the last 
line of HJ are important, in distinction to the double- well 
case (So, a\ are then the annihilation operators for parti- 
cles in the left and right wells, respectively), where they 
both are exponentially suppressed. In particular, we will 
see that the pair-exchange term cx A3 decides upon the 
class of many-body solutions obtained. 

Let |Z) = \N—l, I) be a TV-particle two-mode state, with 
N — I in the ground and I particles in the excited state. 
Then, all many-particle states can be expressed by = 
Si^o^lO- The energy reads with the Hamiltonian (TJ, 
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where the diagonal coefficient ci 
±At(N - 1){N -I - I) + ±A 2 l(l 
di 



e {N - I) + til + 
1) + \M{N - 1)1, 



+ 2)(J + 1)(JV -l-l)(N- I), while oji = 



y/(N - I) (I + 1). Differentiating the functional (5) with 
respect to ipi, the linear system to be solved reads 

- — {uiipi+i + 0Ji-itpi-i) +ciipi. (3) 

There are three important cases, A3 < 0, A3 > 0, A3 = 0, 
for a given set of interaction couplings Ai,A 2 ,A 4 (which 
enter q), assuming SI > 0. 

To get first insight, and introduce some notions, we 
begin by discussing the simplest case of no pair-exchange 
between orbitals and no coupling of levels, A3 = SI = 0. 
For this configuration, exactly one state |Z) will obtain, 
i.e., a definite occupation of both the ground and ex- 
cited single-particle state, because the Hamiltonian is di- 
agonal in I. The necessary condition for such a frag- 
mented state is convexity of the energy parabola in I, 
%£>r = A\ + A 2 — A4 > 0. The energy is minimized 

IOI tmm — 2 ^ 2(Ai+A 2 -A 4 ) ■ UI1 olIlra 

hand, if the energy parabola is concave, we get a con- 
densate, i.e. all particles accumulating in one of the two- 
mode states |7V, 0) (/ = 0), |0, AT) (I = N), depending 
on which of the two states has lower energy. There is 
however one exceptional case. If the energies of the two 
states are identical, leading, by cq = cjy, to the condition 



N=100, Ai=A 2 = l 




FIG. 1: Numerical solutions of (j4} for the indicated parame- 
ter set, for Az = — 1 (left) and A3 = 1 (right), centered around 
I — N/2 (the Ai are given in an arbitrary overall unit of en- 
ergy). The red squares indicate the even I sector (distribution 
of the (pi), and the blue circles the odd I sector (distribution 
of the $;), using the ground-state superposition \(p) — \<k) for 
clarity of representation. 



61 — eo = {N — 1)(^2 — Ai)/2, concavity leads to the co- 
herent superposition of the macroscopically distinct and 
degenerate many-body states \N, 0) and |0, AT) (1, 0, EH] , 
characterized by number fluctuations of, say, the lower 
single-particle state, AN§ = (N§) - (N Q ) 2 = 0(N 2 ). 

In the case of a still vanishing SI = 0, but with a gen- 
eral set of {Ai}, Eq. ([3]) decomposes into two independent 
equations connecting the even and odd I sectors of ipi , 

Elp2k = C 2 feV'2fe + d 2 k1p2k+2 + d 2 k-2lp2k-2, 
Etp2k+1 = C2k+li>2k+l + ^2/5+1^2/5+3 + d2k-1^2k-l, (4) 

Vfe € {0, . . . , N/2}. We have numerically established 
that (to faster than exponential accuracy in N, i.e. 
cx exp[— TV"] with a > 1) the two resulting ground states 
in the sectors of even and odd I, \<p) = J^i't'i |2Z) and 
1$) = J2i ®l 1 2Z + 1) are degenerate when the continuum 
limit (see below) is valid. We display typical results for 
ipi distributions in Fig.[lJ where we put eo = ei = here 
and in the plots to follow (the ej influencing the location 
of the maximum of the ipi distribution, see Eq. J7]) and 
the discussion following it). 

When the pair-exchange coupling A 3 is positive (Fig.[lJ 
right panel), it causes alternating signs, sign (ipi^i+2) = 
— 1, induced by its occurrence in the energy functional 
(2): The energy is minimized by choosing A31P1 ipi + 2 < 0. 
Using the same argument, A3 negative (Fig.[lJ left panel) 
requires the signs of the V/ to be identical. The energy 
contribution of the pair-exchange term is generally neg- 
ative, and its maximization forces ip% to be real. 

The existence of sign changes has profound effects on 
the properties of superposed ground states. The first 
order coherence g\ = ^(a' ai + a\ao) vanishes for any 
superposition a \ <p)+b 1$) in the presence of oscillations in 
the even and odd sectors of ipi, i.e. for A3 > 0, leading to 
well-defined fragmented ground states, in the sense that 
their degree of fragmentation $ in |(5|) is independent of 
the concrete realisation via the weights a and b. To the 
best of our knowledge, this mechanism of fragmentation 
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for positive pair-exchange coupling in a trapped scalar 
Bose gas has not been discussed before in the literature. 
For A 3 negative, on the other hand, the possibility of 
superpositions allowed by the ground state degeneracy 
causes the system to be arbitrarily tunable between fully 
fragmented and coherent states; this can only be avoided 
by lifting the degeneracy by a small but finite Q. 

We define the degree of fragmentation to be £ = 
1 — |Ao — Xi\/N, where Ao,i are the eigenvalues of the 
single-particle density matrix p^) — (aj^a^), given by 

Ao,i/JV = | ± ^-^(JVoiVi-K^P). They col- 
lapse into a single nonvanishing one, A = N, as soon 
as the state becomes coherent, (%ai) = y/NoNi, and 
therefore # = for coherent states. The degree of frag- 
mentation evaluates to 



N=100, Ai=l, A 2 =0.5, A 4 = l, B=l/N 
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where (ajar) = ^_^l(N - I + 1). 

The degeneracy of the two many-body states is lifted 
by fi, which connects ipi and ipi±i, and resolves this de- 
generacy already for Q = 0(1/ N) [so that the energy 
contribution of the fl term is still suppressed to 0(1/ N 2 ) 
compared to the A3 term]. For A3 < 0, turning on fl 
then establishes a coherent state characterized by num- 
ber fluctuations ANq oc N (with a proportionality factor 
of order unity), while for A3 > 0, fragmentation persists 
for £1 > 0. We illustrate the change of # with A3, for 
finite fl, in Fig.[2j We have checked the robustness of the 
fragmentation obtained on the A3 > side upon increas- 
ing the interlevel coupling up to ft ~ O(N ), as well as 
for variations of the values of the single-particle energies 
eo and e\ to the same order. The degree of fragmentation 
is therefore stable for small perturbations on the single- 
particle level, different to what was found in (lil . ItH |. 
where the origin of fragmentation is distinct from our 
interaction-couplings based mechanism. Finally, because 
the energy contribution of the A3 term is negative, for a 
concrete realization with sufficiently small A3 the frag- 
mented state has to be at a (local) minimum of the energy 
in the parameter space of the orbitals ^>i. The latter may 
be determined by a variational ansatz, e.g., by using the 
ellipsoid half-axes of harmonic oscillator trial wave func- 
tions as the orbital parameters 2(J 21 1- 



We now consider the continuum limit of slowly vary- 
ing ipi [4], which allows for an analytical solution in the 
fi ^ limit. Employing the definition j = I — 
and expanding Cj, dj to quadratic order in j/N, neglect- 
ing 0(1/N) terms, we obtain a differential equation for 

\m\, 
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FIG. 2: The degree of fragmentation 5 in Eq. J5]) across the 
transition from coherent to fragmented states at A3 — 0, for 
small but finite f2 = 1/N. The blue dots represent the frag- 
mentation from the numerical solution of (J3}; the solid red 
line is the continuum approximation from in the — > 
limit [the agreement, though not visible on the scale of the 
figure, has also been verified for the A3 < side]. The inset 
shows the first-ord er coherence, with the continuum result 
gi = ^/N 2 /4- 6 2 exp[-l/(4iVv / ^)]. 



Here, E' differs from E by a constant, and the location 
of the distribution center, the shift, is given by 
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N(A 1 - A 2 )/2 + e Q -e 1 
A 1 +A 2 + 2\A 3 \-A 4 ' 



(7) 



The solution of lowest energy of the equation for in 
([6j is a Gaussian centered at 6, cf. Fig.[H and of width 

cr osc = VNR 1/A {{Ai\), provided that 



R{{Ai}) 



I A 3 



A 1 + A 2 + 2\A 3 \ -Aa 



(8) 



is positive. Note that particle number [to 0(1/ N)] and 
single-particle energies do not enter R, in distinction to 
the corresponding criterion in the double- well case 
For either fragmented (A3 > 0) or coherent states (A3 < 
0), a necessary condition besides R being positive, is that 
the shift (O fulfills \&\ <§C y; in particular, a maximally 
fragmented state has 6 = 0. In the continuum limit, 
using the Gaussian shape of the ip(j) distribution, it may 
be shown that the degree of fragmentation is given by 
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A3 > 0, 



1 — e tN~7S\ 1 



(iV/2) 



1 — e A3 < 



(9) 



The above analytical approximation accurately repro- 
duces the numerically obtained result, cf. Fig. Note, 
provided we have the scaling e\ — eo = CN 1 for large TV, 
with 7 < 1 and C constants, assuming that all Ai scale 
identically and A\ ^ A 2 when C ^ 0, we conclude from 
Eq. ([7]) that fragmentation persists for N — > 00. This is 
a remarkable feature of the present two-mode model. 

In case that R becomes negative, or diverges, in the 
nondegenerate case Co 7^ cat we simply obtain Fock states 
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corresponding to the lower two-mode state, i.e. either 
\N,0) or |0, N). The degenerate case leads to coher- 
ent superpositions of macroscopically distinct quantum 
states, for which the original single peak in the ipi distri- 
bution splits into two located around I = and I = N. 
The crossover is illustrated in Fig.[3l by varying A4 in 
((H|, approaching the divergence of R. It appears that 
a fine-tuning of parameters is necessary, illustrating the 
sensitivity of the superposed states to parameter fluctu- 
ations [151 ]. 

N=100, A 1 =A 2 =A 3 = l, A 4 e [3 . 92, 4 . 02] 
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FIG. 3: Crossover to a coherent superposition of degener- 
ate macroscopically distinct quantum states for diverging R 
[Eq. iJH}], using a maximally fragmented state with © = 0, 
upon variation of A4, across the pole of R. Dark blue to light 
green (dark to light gray) indicates increasing A4; = 0. 

In conclusion, we have shown that adjusting the four 
interaction couplings in a two-mode model by employing 
the inhomogenity of the single-particle orbitals in a trap, 
in conjunction with a generally anisotropic and nonlo- 
cal interaction, leads to fragmented ground states whose 
degree of fragmentation is not sensitive to perturbations 
on the single-particle level. The physical origin of this 
interaction-based fragmentation is the pair-exchange pro- 
cess of bosons in a trap, not occurring in the continuum. 
Pair exchanges cause oscillations in the distribution of 
the many-body wavefunction amplitudes, which as a con- 
sequence entail vanishing first-order phase coherence for 
positive pair-exchange coupling. 
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grant No. FI 690/3-1. 



[1] O. Penrose and L. Onsager, Phys. Rev. 104, 576 (1956). 
[2] C.J. Pethick and L. P. Pitaevskil, Phys. Rev. A 62, 
033609 (2000). 

[3] P. Nozieres and D. Saint James, J. Physique 43, 1133 
(1982). 

[4] R.W. Spekkens and J. E. Sipe, Phys. Rev. A 59, 3868 
(1999). 

[5] D. Jaksch et al, Phys. Rev. Lett. 81, 3108 (1998); M. 

Greiner et al, Nature 415, 39 (2002). 
[6] For a summarizing review, see E. J. Mueller, T.-L. Ho, M. 

Ueda, and G. Baym, Phys. Rev. A 74, 033612 (2006). 



T.-L. Ho and S. K. Yip, Phys. Rev. Lett. 84, 4031 (2000). 
J. Dukelsky and P. Schuck, Phys. Rev. Lett. 86, 4207 
(2001). 

O. E. Alon and L. S. Cederbaum, Phys. Rev. Lett. 95, 
140402 (2005); L. S. Cederbaum, A.I. Streltsov, and 
O. E. Alon, Phys. Rev. Lett. 100, 040402 (2008). 
N. K. Wilkin, J. M. F. Gunn, and R. A. Smith, Phys. Rev. 
Lett. 80, 2265 (1998). 

See also the discussion in [2| of the example for fragmen- 
tation given in [Toj |. 

G. J. Milburn, J. Corney, E. M. Wright, and D. F. Walls, 
Phys. Rev. A 55, 4318 (1997). 

Such superposed macroscopic quantum states are some- 
times referred to as "Schrodinger cats", e.g. in For 
clarity, we retain the explicit description of these states. 
T.-L. Ho and C.V. Ciobanu, J. Low Temp. Phys. 135, 
257 (2004). 

J. I. Cirac, M. Lewenstein, K. M0lmer, and P. Zoller, 
Phys. Rev. A 57, 1208 (1998). 

R. Kanamoto, H. Saito, and M. Ueda, Phys. Rev. Lett. 
94, 090404 (2005). 

A.D. Jackson, G. M. Kavoulakis, and M. Magiropoulos, 
Phys. Rev. A 78, 063623 (2008). 
A.J. Leggett, Rev. Mod. Phys. 73, 307 (2001). 
A possible realization of finite SI is to use two hyper- 
fine states coupled by a two-photon Raman transition. 
Counterpropagating laser beams with effectively space- 
dependent Rabi frequency fiR a bi(r) lead to different spa- 
tial orbitals for the effective internal states [C. Zimmer- 
mann, private communication]. 
[20] For an illustrative example, we use as the two states 
the ground and an excited state of a harmonic oscillator 



[10 

[11 
[12 
[13 

[14 

[15 

[16 

[17; 

[18 
[19 



-p 2 /2a 2 

„2 /o„2 



z 2 /2b 2 ] 

-2 / 7,21 



along z, ^o(p,z) = (aVbiri) 1 exp [- 

and *i(p,z) = {ab^/b/^Tvi'^zexp [-p 2 /2a 2 - z 2 /2b 2 
(assuming cylindrical symmetry, where a and b are the 
half-axes of an ellipsoidal cloud [21]). Writing the dipole- 
dipole interaction as V(r) = 4 ^, 3 (1 — 3« 2 /|r| 2 ), where 
gd is the dipole coupling constant [22], and adding the 
usual contact term, we obtain A\ = ^ V 3 \g + ffj/i (cr)] , 

A2 = §(^[S + 9*M<0]. A 3 = l T7 ^[ g + g dfs (a)], 
and A4 — 2A\, where the /»(<x) are functions de- 
pending on the aspect ratio a = a/6. They allow to 
sample the anisotropy and long-range character of the 
dipole-dipole interaction over the whole range from 
oblate (a > 1) to prolate (a < 1), while with contact 
interactions only, the ratios of the Ai are all constant 
Hi]. Explicitly, /1 = [2cr 2 + l-3a 2 ath(cr)], f 2 = 

[2a 6 + 4a 4 + 8cr 2 + 1 - cr 2 (7 + 4a 2 + 4cr 4 )ath(o-)] 



and / 3 



(<T 2 -1) 2 



[2cr 4 - 10a 2 - 1 + 9a 2 ath(a)] , where 



ath(a) = arctanb/c 2 — 1]/Vc 2 — 1- 
[21] 0. Elgar0y and C. J. Pethick, Phys. Rev. A 59, 1711 
(1999). 

[22] U.R. Fischer, Phys. Rev. A 73, 031602(R) (2006). 

[23] Going beyond the dominance of contact interactions in 
monatomic gases, by tuning the relative contributions 
of contact and dipole-dipole interactions using Feshbach 
resonances, has been experimentally demonstrated in T. 
Lahaye et al, Nature 448, 672 (2007). 

[24] We observe here that the authors of [2ll ] were in fact 
considering what we call shifted coherent states, which 
have only one eigenvalue of the single-particle density 
matrix, instead of fragmented states. 



